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ABSTRACT 



During periods of increased magnetospheric activ- 
ity, whistler wave emissions have been observed with 
increased wave amplitudes. We consider a pitch an- 
gle diffusion mechanism that is shown to scale with 
whistler wave amplitude and hence is 'switched on' 
during these periods of intense activity. 

We consider the interaction between relativistic elec- 
trons and two oppositely directed, parallel propa- 
gating whistler mode waves. We show that for in- 
tense whistlers this process is stochastic and results 
in strong electron pitch angle diffusion. 

We show that the interaction is rapid, occur on 
timescales of the order of tens of electron gyrope- 
riods and that the interaction is sensitive to wave 
amplitude, wave frequency and electron energy. 

Key words: Electron, Whistler, Substorms, Diffu- 
sion. 



1. INTRODUCTION 



Electron whistler interactions have long been consid- 
ered as a mechanism for pitch angle scattering. Gy- 
roresonance processes with near-parallel p ropagating 
waves (e.g. [Kennel and Pctschek (1966| ), Lyons anc 
Williams (1984| ) j have been shown to produce pitch 
angle diffusion for electrons that are at resonan ce 
with a series of whistler waves (see Gcndrin (1981 )). 
We generalize resonant diffusion to include all phase 
space dynamics, i.e. as well as considering the reso- 
nant diffusion of trapped electrons we consider diffu- 
sion of untrapped electrons (we refer to this process 
as 'off- resonance' diffusion). Therefore we maximize 
the area in phase space contributing to pitch angle 
diffusion. 



The underlying dynamics of the interaction between 
electrons and a single w histler are in herently sim- 
ple, as demonstrated by Laird (1972 ) who derived 
the Hamiltonian for relativistic electrons interacting 
with a whistler with a single wave number. However, 



for a single monochromatic whistler the process is not 
stochastic. We introduce stochasticity by including 
an identical, oppositely directed whistler mode wave. 

We initially consider a simplified system consisting 
of monochromatic whistlers in order to understand 
the underlying behaviour. This treatment is then 
extended to consider whistler wave packets, i.e. a 
group of waves with a range of frequencies. 

We derive approximate equations in the limit of low 
whistler wave amplitudes and consider the degree of 
pitch angle diffusion for waves of different frequencies 
and bandwidths and for electrons of different ener- 
gies. 



2. EQUATIONS OF MOTION 



We derive full relativistic equations of motion and 
approximate them in the limit of low wave ampli- 
tudes, for the case of monochromatic whistlers and 
wave packets. We normalize time to the electron 
gyrofrequency fl e = eBo/m e , (where Bq is the back- 
ground magnetic field), the wave amplitude is nor- 
malized to the background magnetic field, wave fre- 
quency is normalized to the gyrofrequency and we 
normalize the electron velocity to the phase velocity 
of the waves, given by the electron dispersion relation 
(ignoring ion effects): 
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where u pe is the plasma oscillation frequency. Elec- 
trons can become trapped by either of the two 
whistlers if they are traveling close to the resonance 
velocity, given by: 



k ■ V r 
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where n is an integer and 7 = (1 — v 2 /c 2 ) 1 ' 2 is the 
relativistic factor. 
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Figure 1. Stroboscopic surface of section plots for the monochromatic whistler interaction. For quiet time wave 
amplitudes, panel a), all trajectories are regular, with resonances given by the resonance condition Equation 
For high wave amplitudes observed during disturbed times, panels b), phase space is dominated by stochastic 
trajectories with regular trajectories confined to close to the resonances. The stochastic region is bounded above 
and below by the first regular, untrapped, trajectories away from resonance, therefore there is a limit on the 
fusion of electrons in phase space. 



2.1. Monochromatic Whistlers 



In the limit of low wave amplitudes the full equations 
of motion can be reduced to: 
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where x is the distance along the background mag- 
netic field, b is the normalized wave amplitude, 
vj_q is the initial perpendicular velocity and 70 = 

lMl-uxg/c 2 ). 



fx 



2&w_lo . 
= \ l e sin 

7o 

. n 2£w±o . 
+ \ l e sm 

7o 



( uj)t+kox 

.70 

( u>)t — kox 

.7o 



sm[(t-[3x)Auj/2} 



(t - fix) 
sin[(t+/3a;)Aw/2] 



(t + (5x) 



(5) 



where 1//3 = duj/dk is the group velocity of the 
waves. The wave packet equation (^) yields the 
monochromatic whistler equation (^|) in the limit 
Auj — > 0, with amplitude b' — btt e Auj. 



3. NUMERICAL RESULTS 



2.2. Wave Packet Approximation 



Instead of a pair of waves it is more realistic to con- 
sider the interaction of a wave packet, ie a small 
group of waves with a range of wave frequencies and 
wave numbers. We assume the wave amplitude is 
non-zero over wave frequency range Auj: 



f ° 

b(ui) = { b 
{ 



lo < luq — Alu/2 

ui - Au/2 < uj < loq + Aiv/2 (4) 
u) > LJ{) + Auj 



where luq is the central wave frequency of the 
wave packet. Integrating the monochromatic ehistlcr 
equation over the frequency range Auj gives the 
following wave packet equation: 



The monochromatic and wave packet equations were 
solved numerically using a variable order, variable 
stepsize differential equation integrator. We con- 
sider physical parameters for the terrestrial magne- 
tosphere at L — 6: gyrofrequency, Q e — 25.3 kHz, 
plasma frequency, uj pe = 184 kHz, background mag- 
netic field, Bo = 144 nT and wave amplitude = 
0.5 pT, giving a normalized wave amplitude, con- 
sistent with quiet times in the terres trial magneto- 
sphere , _ b = 10~ 5 ( see f or example |Nagano et al. 



(1996D , [Parrot (1994| ) and [Bummers and Ma (200q )^ 



The phase plots in Figure |l] are comprised of stro - 
boscopic surfaces of section ( |Benettin et al., 1976 ) 
to sample the full electron phase space. The initial 
parallel velocity was varied over the range [—v r , v r ], 
where v r is the resonance velocity, given by the res- 
onance condition, (@), for n — 1. All electrons were 
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given a constant initial perpendicular velocity, with 
u_L rs 20v r as it was found that a high velocity 
anisotropy was required for stochasticity. 

In Figure [l] we plot parallel velocity uy = dx/dt 
against phase angle ^, where ?/? is the angle between 
whistler propagating in a positive direction along 
the background field and the perpendicular veloc- 
ity, v±. In panel a) we consider a whistler wave 
amplitude consistent with quiet times (b = 1CP 5 ) 
all trajectoriesin phase space are regular. There is 
little change in vn and hence only weak pitch an- 
gle diffusion. As the wave amplitude is increased, 
stochastic trajectories are introduced, as the regular 
trajectories between the two resonances are progres- 
sively broken down. In panel b) we consider the case 
of intense whistler wave activity durin g su bstorms 
(b = 10~ 3 , se e for example Parrot (1994) and Nagano 
et al. (1996)). The stochastic region grows to en 



so 



compass the resonances as the wave amplitude is in- 
creased. Regular trajectories are confined to KAM 
surfaces (near-integrable traject ories with an approx- 
imate constant of the motion ( Tabor, 1989 )). The 
stochastic region is bounded by the first untrapped 
(regular) trajectories away from the resonances, thus 
there is a limit on diffusion in phase space. 

As well as resonant diffusion of trapped electrons, 
there is diffusion of untrapped electrons throughout 
the stochastic region of phase space. Since, for suffi- 
cient wave amplitudes, the stochastic region can en- 
compass the resonances, the diffusion of untrapped 
electrons, which we refer to as 'off-resonance' diffu- 
sion, may be enhanced over resonant diffusion. In ad- 
dition we achieve pitch angle diffusion from a larger 
area of phase space. 

Due to the time dependent nature of the wave packet 
equation (||) it is not possible to construct phase dia- 
grams as in Figure ^ for the monochromatic whistler 
case. Instead we can consider the dynamics of single 
electrons. In Figure ||we show a single trajectory so- 
lution of the wave packet equation rtq) , for quiet time 
wave amplitudes and wide (Auj = fi e /50 = 500 Hz) 
and narrow (Auj = Sl e /500 = 50 Hz) whistler 



wave packets (see for example Carpenter and Sulic 
(19p8[)). We consider the change in pitch angle from 
an initial pitch angle of 90°. For narrow wave packets 
there is little change in pitch angle and the trajec- 
tory is regular. For wide wave packets the trajec- 
tory is stochastic with a large change in pitch angle 
(Aa ~ 25°) occuring within a few tens of electron 
gyroperiods (~ 10 ms). We can now achieve strong 
pitch angle diffusion for wave amplitudes consistent 
with the quiet time magnetosphere. 



4. PITCH ANGLE SCATTERING 



Using the wave packet equation (||) with quiet time 
whistlers (b = 10~ 5 ) and a relati vely wide band 
whistler (A uj — f2 e /50 = 500 Hz, ( Carpenter anoj 
[Sulic, 198S} )), we can estimate the degree of pitch 
angle scattering. In Figures and we estimate the 
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Figure 2. Change in pitch angle from an initial 
pitch angle of 90°, for quiet time wave amplitudes 
(b = 10~ 5 ) and narrow and wide whistler wave pack- 
ets. For narrow wave packets there is little change in 
pitch angle. For wide wave packets there is a large 
change in pitch angle (Aa ~ 25°) occuring within a 
few tens of electron gyroperiods (~ 10 ms), hence the 
interaction is rapid. Changes in pitch angle attenu- 
ate with time and the pitch angle reaches a constant 
value. 
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Figure 3. Log change in pitch angle, (log 10 |l+Aa|), 
as a function of wave frequency, u>, and initial pitch 
angle, oto, for 10 keV electrons (v = 0.2c) and quiet 
time (low amplitude) , wide band whistler wave packets 
(Auj = Q e /50). For high to moderate initial pitch 
angles (oiq = 50° — 90°^ there is a change in pitch 
angle of up to 40°. For low pitch angles (ao = 5° — 
10° j the change in pitch angle is of the order of a 
few degrees. 
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log change in pitch angle, (log 10 |1 + Aa|), as a func- 
tion of wave frequency, u>, and initial pitch angle, ao- 
We consider the interaction between lOkeV electrons 
in Figure |H (lOOkeV electrons in Figure ||) and wide 
band whistlers (Auj — ui e /50). For high to moder- 
ate initial pitch angles (ao = 50° — 90°) there is a 
change in pitch angle of up to 40°. For low pitch 
angles (ao = 5° — 10°) the change in pitch angle is of 
the order of a few degrees. In Figure |] we see a sim- 
ilar degree of diffusion except that lower frequency 
whistler wave packets are required. 

Change in pitch anglefci 100 keVeections 




Initial Pitch Angle, d Q 

Figure 4- As for Figure |^ except for 100 keV elec- 
trons (v = 0.5c). Pitch angle scattering from a few 
degrees up to 40° occurs although lower freguency 
whistler wave packets are required. 



5. DISCUSSION 



We have considered electron-whistler wave particle 
interactions to investigate diffusion over all phase 
space, to include both resonant and 'off-resonance' 
diffusion. We have considered a simplified interac- 
tion with monochromatic whistler wave to under- 
stand the underlying behaviour and have shown that 
the presence of the second whistler wave introduces 
stochastic effects into the system. For wave ampli- 
tudes consistent with disturbed times we have shown 
that 'off-resonance' diffusion occurs and that reso- 
nant diffusion is unchanged. 

We have considered a more realistic case of whistler 
wave packets and have shown that for relatively wide 
band whistler wave packets strong pitch angle dif- 
fusion occurs for wave amplitudes consistent with 
quiet, undisturbed, times. For high initial pitch an- 
gles we estimate a change in pitch angle of up to 40°, 
while for low pitch angles a change of a few degrees 
is estimated. 

The effectiveness in scattering electrons of different 
energies is dependent on the wave frequency. Elec- 
trons with low energies (10 keV) are readily scattered 
by waves of around half the electron gyrofrcquency, 



while electrons at higher energies (100 keV) are scat- 
tered by lower frequency wave (cj ~ f2 e /10). MeV 
electrons would require extremely low frequency 
waves for efficient scattering, hence our mechanism 
is most efficient for electrons in the 10 — 100 keV 
range. 
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